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Abstract 

We discuss 1/Q corrections to hard processes in QCD where Q is a large mass pa- 
rameter like the total energy in e + e~ annihilation. The main problem we address 
ourselves to is whether these corrections to different processes (concentrating for 
definiteness on the Thrust and the Drell-Yan cross section) can be related to 
each other in a reliable way so that the phenomenology of the 1/Q corrections 
can be developed. We derive first the relation valid to lowest order using both the 
renormalon and finite-gauge-boson mass techniques to check its independence on 
the infrared cut- off procedure. We then argue that the 1/Q corrections are due 
to soft gluons which factorize into a universal factor such that the lowest order 
relations are preserved in higher orders. 
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1. Power like corrections, (Aqcd/Q) 5 , where Q is a large mass scale provide important 
insight into non-perturbative dynamics in QCD. The best known example are the QCD sum 
rules [0] where the leading power corrections are proportional to the gluonic condensate, 
(G 2 )/Q 4 . The knowledge of the power corrections in that case is based on the OPE. Since 
OPE is essentially a Euclidean construction the limitation of this approach is that it is ap- 
plicable only for physical quantities that can be obtained by analytical continuation from 
Euclidean space, such as the total cross section for e + e~ — > hadrons. The power like correc- 
tions can also be studied H by means of infrared renormalons since the latter are particular 
graphs sensitive to large distances. Thus, all the general properties are observed and the 
nonperturbative contribution results in an overall factor which can be large numerically. 

The advantage of the renormalon technique is that it can be applied to any infrared 
safe quantity and it is by using this technique and its variations that 1/Q corrections were 
discovered recently |3|, f|, ^. |^. |7|, Q . One of the modifications of the renormalon technique is 
introduction of finite gluon mass (A) in the evaluation of the lowest order radiative corrections 
||, 0, [7J. In this case one looks for aX/Q corrections and, in particular, in reference J7J an 
attempt was made to relate the 1/Q corrections to different event shape variables. 

2. Since 1/Q corrections are the leading ones in most cases, their theory is of immediate 
interest. Moreover, the very fact of existence of such corrections follows from simple dimen- 
sional estimates of the infrared contribution to one or another inclusive quantity by means of 
perturbative graphs. The non-perturbative effects bring then an unknown overall factor. To 
develop a phenomenology of 1/Q corrections one needs however to derive relations between 
these corrections to different processes in a reliable way. 

In the absence of the OPE there are apparent difficulties in deriving such relations. In 
particular in the method of introducing a non-zero A it is not at all clear how to go beyond the 
leading order because of the non-Abelian nature of the interaction of the gluons. Similarly, 
in the renormalon technique, (here we follow the approach of ||) higher order corrections 
are unsuppressed because the coupling constant appears to be normalized at a low scale 
(see later discussion). Also, it is far from being obvious that relations between the 1/Q 
corrections obtained in one scheme of introducing an infrared cutoff will be valid universally. 

These are the questions to which we address ourselves in this paper. For sake of def- 
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initeness we concentrate first on a relation between 1/Q corrections to Thrust and to the 
Drell-Yan process. It is straightforward to derive such a relation in lowest order in a and we 
check its validity by deriving it in the two schemes mentioned above. Next we discuss the 
origin of the 1/Q corrections and argue for their universality. The key observation is that 
the same kinematical region of soft gluons is responsible for 1/Q correctins in Thrust and 
the Drell-Yan cross section. The soft part factorizes and the region of the large coupling 
does not affect the relation between the quantities. It is natural, then, to assume then the 
non-perturbative corrections due to the large distances factorize as well. We conclude with 
discussing the other event shape variables. 

3. We proceed now to the computation of the 1/Q correction, to the lowest nontrivial 
order, to the Thrust variable. This quantity is defined as: 

rri ^i Pi ,n\ /-, \ 

T = max n — — - (1) 
^ilPil 

where, pi >n is the component of i-th particle momentum, in the reference direction n which 
is chosen to maximize the right hand side of above. Thus, for a twojet event, T = 1 and for 
a completely symmetric 3-jet event T = 2/3. 

For e + e~ — > qq, T = 1 and in the first order, i.e. for qqg final state 

T = max(xi,x 2 ,x 3 ) (2) 

with Xi = 2Ei/Q in the CM frame where the momenta of q, q and g respectively are denoted 
by pi,P2,P3 and Q is the total CM energy. In the invariant basis 

Q* x . = Q2_ ( g_ ft)2> (3) 

For massless quarks and gluons one finds the range of the variables as: 

< X2 < 1; 1 — x 2 < x\ < 1 and x\ + x 2 + x$ = 2. (4) 

For a nonzero gluon mass A, the above range is modified to 

n ^ ^ i A i A ^ ^ l-x 2 -X/Q f . 

< x 2 < 1 - —, 1 - x 2 - — < x x < . (5) 

Q Q 1 - x 2 
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The cross section for real emission of qqg can be calculated in a straightforward manner in 
terms of o"o, the Born cross section for e + e~ — ► qq. For A = 0, it is: 

1 d<7 0> s Xa + In / n \ 

7T C FTi \T, T" ( 6 ) 



aodx\dx2 2n (1 — xi)(l — X2) 

The calculation of the 1/Q correction to (1 — T) for nonzero A has been reported in ref 
@, so we merely point out some salient features. From (|}|) one first computes the Thrust 
distribution l/a da/dT by substituting the phase space region according to Xi,x 2 ,x 3 being 
the largest fraction and finding 1/cq • da/dT in each region. Adding these, gives the total 
thrust distribution from which (1 — T) is easily obtained. The result is 

(l-T> 1/e = -4C^A. (7) 

We wish to emphasize that the contribution to the 1/Q term only comes from max (xi,^) 
(i.e. the region of energetic quarks + antiquark only) and, in these regions, only from terms 
linear in A coming from the phase space constraints (|5]). It is therefore seen to be soft-gluon 
dominated, i.e. coming from gluon energies of order A. 

At this point, to set the stage for subsequent developments, it is appropriate to give precise 
definitions of what we call soft and what we call the hard collinear regions throughout this 
paper. If we denote a typical momentum by k, and if by 5 we denote a small parameter whose 
vanishing give infrared sensitivity (For example, 5 could be 1 — z, the longitudinal fraction 
of the gluon energy when we consider the r — > 1 limit in inclusive leptoproduction)then the 
two regions are defined thus. 

Soft : k + ~ k„ ~ k± = 0{Q5). 
Hard collinear : k + = 0(Q5), k± = 0(QV5), k_ = 0(Q). (8) 

In the Hard collinear case we take the original hard particle to be moving in the "— " direction. 

In the second method of evaluating the 1/Q corrections to thrust, it will be convenient 
to introduce the variable " spherocity" , S 0] 

S = min (i) 2 (ML) 2 (9 ) 



it J v s *IpI 
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where, is the transverse momentum with respect to the minimum direction. Thus S = 
for 2-jet event. For a three parton final state 

64 

S = — - )(l-x 2 l-x 3 . (10) 

If by x± = we denote the fractional transverse momentum of either parton in opposite 
hemisphere of the fastest parton, then x\ = ^S. Hence 



^ = ^(l-x 1 )(l-x 2 )(l-x 3 ). (11) 



Furthermore, one has the limits: 



1-T) 2 (2T-1) k 2 l-T 

< 7^ < — i— ■ (12) 



T 2 ~ Q 2 ~ 4 

We next change variable: xi,x 2 — > T,kj_/Q 2 = x and from (H) we find for the double 



differential cross section 



1 da a s T 
-C F - 



a dxdT 2tt ^(1-T)(1-4x/(1-T)) 1 /2 
2(T 2 + A 2 + ) 2(T 2 + A 2 ) (2 - T - A+) 2 + A^ (2-T-A_) 2 + A 2 



V(l-T)(l-A+) + (1-T)(1-A_) + (T + A+-l)(l-A+) + (T + A_l)(l - A_) J (13) 
where, 

A ± = l-kfl±^/l^^l- (14) 



Integrating over x gives the well known expression [11] 

Ida a s ^ ( 4 , 2T - 1 ^ r 2T - 1 3(3T - 2)(2 - T)\ 

= — C F — r Zn h 6/n - . (15) 

a dT 7T F \ri-T l-T l-T l-T J v; 



we are interested in 



1 .„ r 1 ^ , I da 



(l-T) = dT , dx (l-T). (16) 

v 1 J2/3 7(i-t)2(2t-i) aQ dxdT K ! y 1 



To get terms linear in Q, we interchange the integration procedure 



1 r 1 ^- r*2 r T 2 



dT I >2 dT ^ dx dT (17) 



and the new limits x±, x 2 and Ti, T 2 are 

Xl = 0,x 2 = -^- T 1 = l-v^+..., T 2 = l-4z + ... , (18) 

where 7\ and T 2 are given as an expansion in x. Proceeding thus, we find that the only 
contribution to the 1/Q term comes from the lower limit of T integration: i.e. for T — > 1 
and the soft gluon region. 

A contribution to < 1 — T >uq ha the form: 

n r l/12 A T r l-4x 4 1 

d-T) 1/Q = ^C F [ ^( i ^ dT - Jr - s ^. (19) 
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Precisely this term gives the ln ^~^ behaviour in the thrust distribution above. Next, to find 
the position of the leading infrared renormalon we allow a s to run according to a — > a s (k\) 
(see, e.g., ref [O]) and in the above integral expand in kj_/Q. The leading term is 



The coefficient of 4— Cp is the same as for the A 7^ calculation, if we identify A with: 



A 



The leading IR renormalon at l/2&o can be explicitly seen by introducing the represen- 
tation m 

a 8 (A£) = / daexp{-a(b ln-±)) (22) 

above and doing the k± integration to get a pole at 2abo = 1. 

We next turn to the 1/Q correction to the Drell-Yan process which was discussed in 
||. H. Following ||, the quantity of interest is the Drell-Yan cross section cI<jdy(t,Q 2 ) / 'dQ 2 
normalized to the structure functions of deeply inelastic scattering Fdjs(x, Q 2 ). In particular, 
we will be interested in the large N moments of this which is known to factorize to all orders 
g Thus: 



l ^N^^drr N ^ daDY ^ Q2) = a Q R(a s (Q 2 ))E(N,Q 2 )F 2 DIS (N,Q 2 ) (23) 



dQ 



5 



with r = Q 2 j s, Q 2 is the invariant mass of the dilepton pair and R(a s (Q 2 ))is a finite function 
as iV — > oo. In the above, E(N, Q 2 ) is the contribution of soft gluons only. For large N we 
are in the kinematical domain r — > 1 and here, E(N, Q 2 ) has large perturbative corrections 
which in the lowest order are of type log 2 N and logN coming from the region of soft gluons. 
It is well known |13| that these exponentiate. In E(N, Q 2 ), the dominant term at large N in 
the lowest order is (to detect the renormalon we need to have the coupling running so that 
we have let a s — > a s (k±), which is the appropriate scale for a s , [[12|, [H 



^/^i^/'-^^M)^. (24) 

J 1 — Z J(l-z) 2 Q 2 k\_ 71 

By interchanging the and z integrations we can identify the 1/Q behaviour which was 
found in PI 



Q JO Kj_ Jl-k±/Q \ 7T 



% N -Vl"^. k Jc^<m. (25) 



Q J0 kj_ \ 7T 

The coefficient of the term proportional to (N — 1) comes from real emission diagram and is 
exactly the same as the 1/Q corrections in (1 — T) and is again related to the soft behaviour 
as r — ► 1. This in turn can be traced back to the maximum allowed transverse momentum 
in deeply inelastic scattering and in Drell-Yan: 

(*i)£2 = | = (26) 

and 

(*i)£L = = Q \ z z? = C 1 -^!- ( 27 ) 

The 1/Q corrections come from the lower limit of the z integration in the first of equation 
( fZ5D which is governed by the Drell-Yan soft dynamics as discussed above. Further, the 
contributing region is the one where the invariant mass of the gluon radiation is vanishingly 
small, of order (1 — r)s as r — > 1. 

The same 1/Q term can also be obtained by the method of introducing a gluon mass A. 
In this case again the only contribution to the 1/Q term comes from the specification of the 
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integration boundaries. The coefficient of this term, linear in gluon mass, in the corresponding 
expression for E(N, Q 2 ) is identical to the result using the other procedure with the same 
identification of A as for the case of Thrust. Furthur, it is easily checked to be coming from 
the region of soft gluons alone. 

4. A few comments are in order concerning these results. Although we have checked 
the independence of the relation between the 1/Q terms on the choice of infrared parameter 
to lowest order, we should look at the problem in higher orders as well. The point is that 
higher orders are not suppressed since a s occurs normalized at a low scale. To demonstrate 
the lack of suppression of higher orders within the renormalon technique, we notice that in 
the leading log approximation 

= ( a *-2mX^) a - (kl) - (28) 

Since the renormalon contribution is linear in Aqcd we conclude that a;?(Ar[). contribution 
is not suppressed. Thus we are led to investigate the origin of the universality of the 1/Q 
corrections obtained in the lowest order. 

Let us discuss the nature of the higher order corrections to the IR safe, event shape 
variables y that vanish in the limit of the 2-jet process (like 1 — T in particular). From the 
lowest order analysis we have seen that the 1/Q corrections come from the region T — > 1 
and of energetic quarks (antiquarks) and soft gluons. Indeed, we have seen that as far as 
these 1/Q corrections are concerned we have in the CM frame 

i-r»of+o§- (») 

where M 2 ps 2p q ■ hi in the limit of soft gluons. This generalizes straightforwardly to all 
orders: Let us first argue for the dominance of soft gluons in 1 — T in higher orders. In 
the next order, the infrared sensitive contribution will come if the additional gluon is either 
soft or hardcollinear, real or virtual. Since (1 — T) is an IR safe quantity, the hard collinear 
gluon will cancell against the virtual process and the only IR sensitive contributions that 
will survive are those involving soft gluons. Continuing in this manner, it is seen that to 
all orders the infrared sensitive contributions that will survive are those involving the soft 
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gluons alone. In fact, to all orders, the IR sensitive contributions to 1 — T can come only 
from the region T — > 1 , where we can write 



M? M 2 



Here we have considered two hemispheres, Hi and H 2 , divided by a plane perpendicular to 
the thrust axis, defined P\ and P 2 to be the vector sum of the momenta in these and denoted 
by Mi and M 2 the corresponding invariant masses. It is then easy to see that for small 
1 — T the above is true. In principle, both soft and hard collinear gluons contribute towards 
small invariant masses, however the latter are excluded since they would cause a collinear 
divergence in 1 — T. Indeed, it is well known that the collinear divergences cancell in the 
total cross-section when we sum over final degenerate states. The same is true for < T > 
since the states degenerate with the hard collinear gluon have the same value of T. With 
this in mind we may write in the limit of soft gluons with momenta k, 

Mf = £ rn) (31) 

with, m 2 = 2kjPi. The next step in the argument is that the soft gluons factorize in the two 
jet process under consideration and hence to all orders we may write 

< 1 - T > 1/Q = R T (a s (Q 2 ))E soft (Q) (32) 

E so ft{Q) comes from soft gluons alone and can be expressed in the form 

rQ 2 dk 2 



E. 



soft 



(Q) = l/Qjf -^k ±leik (a s (k 2 ± )) (33) 



leik{ a s{k 2 1 )) contains contributions from all orders and can be obtained in perturbation 
theory from the eikonal approximation. It has the perturbative expansion: 

l eik (a s ) = -C F + {-f\c F K + .... (34) 

71 71 Z 

with, 

K = C A (67/18 - tt 2 /6) - 10/9T R N f . (35) 
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Rt is not IR sensitive and can be obtained as an expansion in a s (Q 2 ). 

The above arguments emphasize that the origin of the 1/Q terms lies in the soft gluon 
radiation. The same is true for other shape variables like the C parameter, which are infrared 
safe and which vanish in the two jet limit. In fact all of these quantities, generically labelled 
as y earlier are governed by the same soft dynamics as they vanish Hl4| . Thus we have in 
general: 

< y > 1/Q = R y (a s (Q 2 ))E soft (Q) (36) 

In the above the quantities R y are y dependant, but the soft functions E so f t are universal. 

We next turn to a discussion of Drell-Yan (inclusive lepton pair production) for large r, 
in higher orders of perturbation theory. As discussed earlier, the factorization in the r — > 1 
limit has been established as given by Eq.fl23"|) , with E(N, Q 2 ) , purely a soft function. The 
1/Q corrections to all orders is of the form j8|: 

F2 Um N ^f dTT N-i da DY {r Q 2 ) = {N _ 1)R{as{Q 2 ))E {Q) (37) 

f dis\ N iQ )°o Jo dQ 2 

The equality on the RHS is, of course, only for the 1/Q terms. In the limit r — > 1 , the 
invariant mass of the quark gluon system which also has the form of Eq.([H]) , becomes 
vanishingly small . For the quantity < 1 — T >, we have a similar situation expressed by 
Eq . (|30D , in the limit of vanishing invariant mass. In fact, the large corrections due to the soft 
gluons will be the same, in the sense that the functions E so f t are universal hence the same 
in Eq.(|33|) and Eq. (|37|) . In || j ei k is called the cusp anomalous dimension. The factor of 
iV — 1 in Eq. (|37|) is what is nescessary to make the correspondance between the two equations 
complete at r — > 1 since the invariant mass of the quark gluon system , is proportional to 
1 — t for the Drell-Yan case. Because of the relation between < 1 — T > and the other event 
shape variables, we see that in general, in, Eq.(^) and Eq.(^), the function E so f t will be 
universal, whereas the perturbatively calculable (for large Q 2 ) infrared insensitive quantities, 
Ry(&s(Q 2 )), and R(a s (Q 2 )) may be different. 

Apart from the situation discussed above, the same jeik also determines the Infrared be- 
haviour of the quark form factor and the so called velocity dependant anomalous dimensions 
in Heavy Quark Effective Theory ||15 |. The universality of the anomalous dimensions of the 
effective currents in HQET was first noted in Jl6| . 
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Up until now we have only discussed the situation in perturbation theory. We now take 
perturbative calculations with a s normalized at low momenta ( as in the functions E so f t ) as 
representative of large distance dynamics in QCD, including non perturbative effects. Thus 
it is natural to expect a factorization similar to that derived perturbatively for the full theory 
as well. If so then the relations discussed between the 1/Q terms in inclusive lepton pair 
production and the event shape varibles < y >, will be valid here also. 

In conclusion, we have presented arguments that a single non perturbative parameter 
describes the 1/Q corrections in Drell Yan and certain event shape variables. The same 
nonperturbative parameter makes its appearance also in heavy quark physics. All of these 
quantities are related to this universal parameter through overall normalizations. Thus 
theory of the 1/Q corrections seems to be ripe for the confrontation with experiments. 

We would like to thank the Service de Physique Theorique, CE-Saclay for hospitality. This 
work was supported in part by the US Department of Energy. 

References 

[1] M.A.Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys., B147 (1979) 385. 

[2] V.I. Zakharov, Nucl. Phys. B385 (1992) 385. 

[3] H. Contopanagos and G. Sterman, Nucl. Phys. B419 (1994) 77. 

[4] I.I. Bigi, M.A.Shifman, N.G. Uraltsev, and A.I.Vainshtein, Phys. Rev. D50 (1994) 2234. 
[5] M. Beneke and V.M. Braun, Nucl.Phys. B426 (1994) 301. 
[6] A.V. Manohar and M.B. Wise, Phys. Lett. B344 (1994) 407. 
[7] B.R. Webber, Phys. Lett. B339 (1994) 148. 

[8] G.P Korchemsky and G. Sterman, Nucl. Phys., B437 (1995) 415 . 
[9] M. Beneke, V.M. Braun, and V.I. Zakharov, Phys. Rev. Lett. 74 (1994). 
[10] H. Georgi and M. Machacek, Phys. Rev. Lett. 39 (1977) 1237. 



10 



[11] G. Altarelli, Phys. Rep. 81 (1982) I. 

[12] D. Amati, A. Bassetto, M. Ciafaloni, G. Marchesini, and G. Veneziano, Nucl. Phys. 
B173 (1980) 429. 

[13] G. Sterman, Nucl.Phys. B281 (1987) 310; S. Catani and L. Trentadue, Nucl.Phys. 
B327 (1989) 323, B173 (1991) 183. 

[14] R.K. Ellis, D.A. Ross and A.E. Terrano, Nucl. Phys. B178 (1981) 421. 

[15] G. P. Korchemsky and A. V. Radyushkin, Phys. Lett. B279 (1992) 359. 

[16] A. F. Falk, H. Georgi, B. Grinstein, M. B. Wise, Nucl. Phys. B343 (1990) I. 



11 



